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Thermal dilepton production from hot QCD 


1. Introduction and observables 


The thermal production rate of off-shell photons, subsequently decaying into on-shell lepton- 
antilepton pairs, is a rich source of information concerning the properties of the hot QCD medium 
generated in heavy ion collision experiments. If invariant masses corresponding to known vacuum 
resonances are avoided, the result can be expected to be relatively free from non-thermal back¬ 
ground effects. Conversely, a focus on thermal modifications of prominent resonances, particularly 
quarkonium states, can in itself be turned into a useful probe of medium properties. In the present 
contribution we concentrate on non-resonant production and in particular on the contribution from 
gluons and three (Nf = 3) deconfined light quark flavours. 

In order to establish notation, let us denote by T the temperature; by k = |k| the total mo¬ 
mentum of the dilepton pair with respect to the plasma rest frame; by ko the pair’s energy; and 

by 

M=^kl-B ( 1 . 1 ) 

its invariant mass. To leading order in the electromagnetic fine-structure constant ote and for mass¬ 
less quarks, the differential production rate per volume reads ^ ^ 


dkod^k 


«b(^o) 

37r3M2 


d{M'^ 





Nf 

lmn^{ko,k) . ( 1 . 2 ) 

i=l 


Here Ug is the Bose distribution, Qi the electric charge of a quark of flavour i in units of the electron 
charge, and ImlTj^ stands for the imaginary part of a retarded correlator (i.e. a spectral function) of 
one flavour, evaluated in an equilibrium ensemble at a temperature T. 

Given that we are interested in making contact with lattice simulations, the retarded correlator 
is best interpreted as an analytic continuation of an imaginary-time correlator (cf. e.g. ref. [Q]). 
Denoting the imaginary-time coordinates by 


X = ir,x), K={o)„k), (1.3) 

with G)„ = InTn, n € Z, the analytic continuation reads 

Ugikofk)^ (1-4) 

where the spinor t/r represents one flavour, and a suitable ultraviolet regularization is needed for 
defining the Fourier transform at short separations. It is interesting to also consider the space-like 


correlator 

/■i/r . r . . 



dte'^'Ty d^x^ {J^{X) 7^,(0)) , X = (x^,z) . 

(1.5) 

Clearly, 


n^{ko,k) = \ dze^’^^V%{(0„,z)\ 

(1.6) 


The last equation suggests that spatial correlations measured with / 0 have a principal relation 
to the real-time dilepton production rate captured by Imir^. 
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2. Different regimes and previous work 


If eq. ( [T!^ ) is addressed within the weak-coupling expansion, the method to be used depends 
on the parametric magnitudes of k and M. In the following we recall the main cases, denoting by 
g = the gauge coupling. For reference, the leading-order result reads 


— ImllR = 


NJM^ 

Ink 


In ■ 


cos 


h(|f 


) 


cos 


I'd-) 


kt- = 


k^Aik 


(2.1) 


The simplest case to discuss is ^ = 0, because then the result only depends on a single kine¬ 
matic variable, M = k^. If M>7iT, then the NLO correction to eq. o is infrared finite and 
small 1^ 1^ ^ (these results have recently been extended to a finite quark mass 0). However, 
the NLO correction increases rapidly as M decreases; for M ~ gT, the correction is of 0(1) and 
needs to be summed to all orders, yielding a large enhancement |^. More recently, it has been 
realized that the original (“HTL”) resummation is not sufficient for capturing all relevant effects 
for M <gT. The correct infrared behaviour, including a transport peak in ImllR/ko whose width is 
^ g^T/ k^ and whose height determines the electric conductivity ^ UgT jg'^, has been worked out 
in numerical form in ref. 0. 

The phenomenologically perhaps most interesting case concerns the production rate of dilep¬ 
ton pairs with a “soft” invariant mass (M ~ gT) but large spatial momentum {k > nT). In this regime 
the NLO-rate has a logarithmic singularity, which is regulated (as indicated below by M — gT) by 
Landau damping of the spacelike quarks mediating t-channel exchange 00: 


-ImllR = ... -h 


a,NcC,T^ 


In 


1 — 2n^{k) 


( 2 . 2 ) 


2 ^\M^gT^ 

where n^ is a Fermi distribution. In addition, there are finite terms which all contribute at the same 
order because of collinear enhancement, and need to be handled through Landau-Pomeranchuk- 


Migdal (LPM) resummation |13] (LPM resummation incorporates HTL resummation in an ap¬ 
proximation valid for k S> gT). In order to avoid double counting, LPM resummation needs to be 


carefully combined with other processes []14|]. 

In a “hard” regime M S> 7iT, Operator Product Expansion (OPE) techniques become applica¬ 


ble [ 15]. The result is available in a closed form up to NLO: 


— ImH^ = 




! + ■ 




ACCaNrCv 


! + ■ 


Ak^ \ k^T 


27^4 


■-hi 




(2.3) 


An y ' An J ' 9 V* ’ J V 

Yet another approach is to carry out lattice simulations. Lattice QCD being formulated in 
imaginary time, with a time coordinate 0<T<l/r,itis however not possible to measure the rate 
directly, but rather a particular transformation thereof: 


pOQ 

GE{r,k)= — ImnR(ko,k) 
Jo K 


cosh(2j — T)ko 


(2.4) 


sinh(^) 

This means that a part of the contribution comes from the “unphysical” domain ko < k. Even though 
it is in principle possible to invert the relation in eq. (^^) for ImnR(ko,k) [Q], in practice large 


systematic uncertainties are induced [16]. A more controlled approach is to insert an analytically 


determined ImnR(ko,k) into eq. (2.4) and compare the resulting GE{n,k) directly with numerical 
measurements (for recent work and references, see ref. [0). 
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3. Methods for recent developments 


We now turn to works whose scope can be summarized as follows: 

(i) The rate ImlTij has been determined up to NLO in a “generic” regime k,M ~ kT, verifying 


the cancellation of infrared divergences and finding in general a small correction [18 


(ii) For a hard momentum k ~ nT there is an additional parametric scale, M ~ across 

which NLO results valid for M> ttT and LPM-resummed results valid for M<gT can be 
interpolated into each other ||^]. 

(iii) It has been suggested that the relation in eq. ( |l.6| ) can be turned into a non-trivial crosscheck 
of dilepton production rate computations, in the sense that the screening masses associated 
with the |z| ^ 1/r behaviour V^y{(0„,z) can be measured non-perturbatively on one hand, 
and computed through an LPM-type resummation on the other. This permits for a direct 
comparison of the two approaches, without analytic continuation [^]. It has also been spec¬ 
ulated, and demonstrated within the AdS/CFT framework, that a direct “analytic continua¬ 
tion” of screening masses (rather than correlation functions) might allow for a determination 


of the electric conductivity [22j. 


(iv) Finally, in an impressive computation, the LPM-resummed analysis of ref. im has been 
extended up to NLO in the soft regime M<gr, verifying the predictions’ stability [^]. 

The methods used in these studies are as follows. For (i), standard (unresummed) perturbation 
theory is sufficient, with the only complication that because of infrared divergences appearing in 
“real” and “virtual” corrections, intermediate stages of the results need to be worked out in the 
presence of an infrared regulator, which cancels in the end [[T8|, [T^. 

For (ii), an essential ingredient is to realize that the result of the soft regime consists of two 
parts, the LPM-resummed one and another part, referred to as 2 —)• 2 scatterings, which also appears 
in the non-resummed NLO expression. So, for a consistent result in the soft regime, the LPM- 
resummed result and the NLO result need to be summed together, but only after subtracting those 
terms from the NLO result which got resummed into the LPM one. This can be summarized as 




expanded to NLO 


-|- imn. 


I full 
rIlpm 


(3.1) 


This expression is correct in the soft regime M<gT because the subtraction removes the danger 
of double-counting terms appearing in Imflpl^pQ, but also in the hard regime M>7rr, because the 
two variants of the LPM expression cancel against each other there (up to higher-order corrections). 

For (iii), two ingredients are needed. One is a theoretical analysis showing that a perturbative 
determination of the screening masses at ft)„ / 0 involves solving inhomogeneous Schrddinger- 
type equations which are just analytic continuations of the corresponding equations appearing in 
LPM resummation (with, in particular, the same “potential” [M, ^ j and the same “inhomogeneous 


terms”). The other ingredient is measuring screening masses with standard lattice QCD; in ref. [ |21[ ] 
this was done for Nf = 2 light dynamical flavours. 

For (iv), fhe paramefers appearing in LPM resummafion, in particular fhe asymptotic fhermal 
masses and fhe pofenfial, need fo be modified fhrough NLO correcfions. Many ingredienfs are 
known from previous work [24, 271, buf need fo be puf fogefher in a consistenf fashion S]. 
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pert. ^-^lat. ^-^pep 

Figure 1: A comparison of screening masses related to different polarization states of the vector current 
in the presence of a Matsubara frequency {(0„ = 2nTn) in dynamical Af = 2 QCD at F ~ 250 MeV (from 
ref. pl|]). The correlation functions were measured in the direction of the z-axis, and were averaged over 
the “transverse” (x^^X 2 ) plane. The good agreement between lattice measurements lat”) and re summed 
perturbative results (“ pert”) suggests that the systematic uncertainties of both approaches are moderate. 


4. Results 


In fig. I], screening masses measured in Af = 2 lattice QCD at T « 250 MeV are compared with 
corresponding predictions from perturbation theory, based on solving an LPM-type Schrddinger 
equation. It is seen that even at this “low” temperature, reachable in the current generation of LHC 
experiments, resummed perturbation theory does reproduce the qualitative features of the lattice 
measurements, with differences only on the ~ 15% level. (It would be interesting to increase the 
resolution of the simulations, to take the continuum limit, and to study several temperatures, in 
order to see if the remaining discrepancies decrease.) 

In fig. ^ a result for ImlTj; based on eq. ( |3.1[ ) is plotted as a function of M and k. A crossover 
from one type of behaviour at M<r to another shape at M>nT is clearly visible; as discussed 
above, the crossover takes parametrically place at M ~ 

Finally, in tig. ^ the physical dilepton rate from eqs. ( |3.1D , ( |1.2[ ) is illustrated. The grey bands 
indicated the uncertainty associated with variations of the renormalization scale. Even though 
the uncertainty is >50% for M< 1 GeV, the overall shape of the curves as well as their general 


normalization can be predicted. We also note that a recent study of the soft regime up to NLO [23] 
increases the results by ^(10%) for general parameters (up to 100% for M = 0.25 GeV), which 
suggests that the grey bands do indeed capture the magnitude of uncertainties. 


5. Outlook 

There are two major “applications” for the results reviewed here. First of all, curves such 
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Figure 2: The vector channel spectral function determined through eq. (3.1) (from ref [Effl). 


T = 0.5 GeV, Nj = 3 



Figure 3: Physical p production rates based on eqs. (|3), (Q (from ref @]). 


as shown in fig. complemented by results for the regime ko < k, can be inserted into eq. 0 
and compared with direct lattice measurements. Compilations suitable for this purpose, based on 
refs. 0,0II, n, 0, can be downloaded from ref. Second, results such as shown 

in fig. 1^ could be inserted into hydrodynamical models such as ref. [29], and compared with 
experimental data from heavy ion collision experiments. Data suitable for this purpose, from 
refs. [pOl p3l], can be downloaded from ref. [pOl]. 
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